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ABSTRACT 

Phase-induced amplitude apodization (PIAA) coronagraphs are a promising tech- 
nology for imaging exoplanets, with the potential to detect Earth-like planets around 
Sun-like stars. A PIAA system nominally consists of a pair of mirrors which reshape 
incident light without attenuation, coupled with one or more apodizers to mitigate 
diffraction effects or provide additional beam-shaping to produce a desired output pro- 
file. We present a set of equations that allow apodizers to be chosen for any given 
pair of mirrors, or conversely mirror shapes chosen for given apodizers, to produce an 
arbitrary amplitude profile at the output of the system. We show how classical PIAA 
systems may be designed by this method, and present the design of a novel 4- mirror 
system with higher throughput than a standard 2-mirror system. We also discuss the 
limitations due to diffraction and the design steps that may be taken to mitigate them. 

Subject headings: methods: analytical, techniques: high angular resolution, instrumen- 
tation: high angular resolution 



1. Introduction 

The challenge of imaging an Earth-like planet around a Sun-like star has two primary compo- 
nents: contrast and angular resolution. An Earth twin orbiting a sun twin at 1 AU at a distance of 
10 parsecs will be 10 10 times fainter than its host star, while being separated by only 100 milliarc- 
seconds. For a 2m telescope with a circular aperture operating at 500nm — a reasonable expectation 
for a visible-band space telescope — this places the planet at 1.95A/-D in the image plane. Here, the 
intensity in the point spread function (PSF) from the host star will be ~ 6 x 10~ 3 lower than the 
peak of the diffraction pattern and over six orders of magnitude higher than the planet, making 



1 Copyright 2012 California Institute of Technology. Government sponsorship acknowledged. 



-2 - 



detection of the planet virtually impossible. To overcome this, a number of optical systems — 
including coronagraphs, occulters, and interferometers — have been proposed to reduce the residual 
stellar light in the regions of the image plane where we wish to look for planets. 

One of the most promising corona graphs for space applications is the phase-induced ampli- 
tude apodization (PIAA) coronagraph (lGuyonll2003l : iTraub k, Vanderbeill2003l : iGuyon et alJl2005l ; 
Vanderbei fc Traubll2005l ) . which uses a pair of mirrors to losslessly remap incident light to apodize 
the beam while retaining high throughput. In addition, incident wavefronts with an angle to the 
optical axis of the system are displaced disproportionately from the optical axis, providing an 
effective magnification of the inner working angle. 



A number of these s ystems have been designed (IPluzhnik et al.ll2006f) a nd tested (jBelikov et al 



2009: 


Kern et al. 


2009: 


Guvon et al. 


2010: 



201ft IPuevo et al 



2011a 



), in 

conjunction with apodizing elements to redu ce manufac t uring requirements on mirrors or mitigate 
diffraction from mirror edges. As shown in IVanderbeil ( 20061 ) . without at minimum an apodizer 
before the mirrors (a pre- apodizer), the effects of diffraction from the edges of the mirrors limit 
the achievable contrast severely, to 10 -5 . Apodizers after the mirrors (a post- apodizer) have also 
proved useful for performing a portion of the beam-shaping, at the expense of some throughput; 
using a post-apodizer can reduce the manufacturi ng requirements on t he mirrors, particularly the 
curvatures required on the edges of some designs (IPluzhnik et al.l 120061 ) . 



While the tools for acc urately modeling propagation through P I AA systems have become fairly 



2006b: 
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Krist et al. 
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Puevo et al. 



2011bl ). the tools for finding new apodizations and mirror shapes are not as developed; in particular, 
the interaction between multiple remapping and apodizing components in the design has not been 
fully explored. In this paper, we discuss the propagation of radially-symmetric amplitude functions 
through PIAA systems, in particular setting down a series of equations in Sec. [2] which allow for 
the arbitrary specification of amplitude profiles throughout the system. This permits mirrors and 
apodizers to be designed explicitly to work jointly as a system; examples of this are shown in Sec. 
[3j Applications to design and modeling of four-mirror systems are discussed in Sec. HI including 
their interactions with apodizers in the system. 



2. Design of amplitude profiles 



A key problem in designing a PIAA system is determini ng what shapes the mirro r s must 
take on to correctly redistribute the light. Early PIAA work (jTraub Vanderbeil 120031 : I Guvon 



20031 ) approached this from the perspective of ray optics, shaping the mirrors to redirect incoming 



rays into a desire d amplitud e distr ibution while maintaining constant pathlength for each ray. 
However, in 2006, 1 Vanderbeil (120061 ) showed that this approach was insufficient; mapping of rays 
was insufficient to capture all of the physics seen in the full diffraction integral. Rays do allow for 
a straightforward visual interpretation of the remapping induced by PIAA, and the workings of a 
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PI A A system using rays are shown in Fig. [TJ (We note that, as an approximation, design is done 
in an on-axis configuration, despite the fact that the system geometry requires that mirrors must 
eventually be placed off-axis; the system layout is generally chosen to minimize the off-axis angle.) 

In the following, we investigate an alternative approach to creating PIAA systems with pre- 
a nd post-apodizers. Rather than ray optics, we start with the Huygens diffraction integral as set out 



m 



Vanderbeil (|2006l ). with a uniform amplitude distribution. We then consider the first term of an 



asymptotic exp ansion of the full H uygens diffraction integral — that is, apply the stationary phase 



approximation (jBorn <fe Wola Il999l ) — and show that the resulting approximation directly relates 
the output amplitude distribution and the mirror shapes, and additionally recovers the geometric 
optics ray-remapping relationship from the full diffraction integral for a uniform input amplitude. 
In Sec. 12.21 we then use the stationary phase approximation to show how to remap non-uniform 
but radially-symmetric input amplitude distributions. 

Our purpose in doing this is threefold: first, to show that a relation between input and output 
amplitude profiles and the mirror heights can be derived from the diffraction integral in a manner 
consistent with, but without recourse to, ray optics, and to be able to make the same sorts of 
statements about the amplitude distributions in a diffractive PIAA system that can be done with 
ray optics. We would note that the purpose here is not to explicitly rederive the ray-optics mapping, 
though the equivalence does introduce a certain amount of convenience. 

Second, the stationary-phase approximation gives a straightforward derivation for the relation- 
ships between the amplitude profiles at different planes that to our knowledge has not been shown 
in a ray-optics context for PIAA previously. 

Third and finally, we wish to lay the groundwork for designs which explicitly minimize diffrac- 
tive effects by appropriate choices of apodizers and mirror shapes. The differences between the 
stationary phase approximation and the full diffraction integral can be expressed explicitly as 
higher-order terms in the same asymptotic expansion; the derivation of these terms falls beyond 
the scope of this current work, but their characteristics and general dependence on system param- 
eters are known and discussed later in Sec. 12.51 along with strategies for mitigating them. 

Sec. 12.31 and 12.41 show how the equations of Sec. 12.21 may be inverted to allow the PIAA 
mirror shapes to be derived directly from specified input and output profiles, as well as how exactly 
the required profile on the input beam is set by the mirror shapes and the output profile. In 
particular, the results of these sections allow one to choose apodizers and a desired output amplitude 
profile — based on whatever criteria which will optimize performance for the specific application — 
and determine exactly what shapes the mirrors must be to make this happen, a novel result for 
radially-symmetric PIAA systems. 
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Fig. 1. — The basic PIAA system, illustrated with rays: light entering the system is redirected 
by the first mirror to concentrate or spread out the light, and the second mirror corrects the path 
length so the wavefront has a uniform phase at the output as well as the input. The mirror depths 
have been increased by a factor of 5 for clarity. Coordinates related to arrangement of these mirrors 
are shown to the right. 
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Fig. 2. — A diagram showing the three amplitude profiles associated with a PIAA system. A m (r) 
is the output amplitude profile when the input amplitude profile to a set of mirrors is uniform, 
while for an arbitrary input profile Ai(r) and the same set of mirrors, the output profile is Af(r). 
As shown later, knowledge of any two of these profiles allows the third to be determined. 
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2.1. The stationary phase approximation applied the diffraction integral 

For a radially-symmetric PIAA system, there are three different amplitude profiles to consider 
when propagating an on-axis wavefront: 

• The amplitude profile of the wavefront incident on the first mirror, which we denote Ai(r). 

• The amplitude profile directly following the second mirror, which we denote Af(r). 

• The 'intrinsic' profile of the mirror pair; this profile would be produced by an incident plane 
wave when considering the remapping as ray-optics, without diffraction. We denote this 
A m {r). 



The amplitude profiles Aj(r), Af(r), and A m (r) are shown in Fig. [5J These three profiles are 
not independent, but rather selecting any two determines the third; in the following we will show 
this and further show how to determine the third given the other two. (A fourth profile from a 
post-apodizer, A p (r), may be included to reshape the amplitude profile further; this factor is solely 
multiplicative and does not factor into the remapping equations that follow. See Sec. 13.11 ) 

Assume we have a PIAA system with intrinsic profile A m (r), assumed real and positive. Con- 
sider initially the output from a wavefront Eq(x, y) incident on a PIAA system with no pre-apodizer; 
following the Huygens-Fresnel principle, the electric field at any point in the output plane can be 
considered to be a superposition of spher ical waves from a ll points across the open aperture, which 
in Cartesian coordinates takes the form (jVanderbeil 120061 ) 

Ea{x,y) 



E(x,y) 



-fQ(x ,y ,x,y) ^ ^ 



(1) 



i\Q(x ,y ,x,y) 

Here, xo and yo are coordinates in the plane of the input aperture, while x and y are coordinates 
in a plane just following the second mirror. In the case of an on-axis plane wave, Eq(x, y) = 1. The 
optical path length Q is defined as 



Q(x ,y ,x,y) = y 7 (x - x) 2 + (y - y) 2 + (z + h(r) - h (r )) 2 + (h(r) - h Q (r Q )) 



(2) 



ri + 



2r$r cos 



h) + {z + h(r) - h (r )) 2 + (h(r) - h (r )), (3) 



where (tq,9q) and (r, 9) are polar coordinates equivalent to (xo,yo), and (x,y), respectively. The 
shapes of the first and second mirror are given by /io( r o) an d h(r), respectively, and defined by 



and 



dh (r ) _ dh(r) _ r (r) - r 
drn dr 2z ' 



ho(0) = h(Q) = 0. 



(4) 



(5) 



We note that this notation has been adjusted slightly from the definitions in IVanderbeil (|2006l ). so 
that both ho and h are between and the maximum height of the mirror. Here, unlike in Eq. [3j 
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?"0 in Eq. [J] is not an independent variable from r; rather, we recall the equations for a single pair 
of radially-sym metric PIAA mirrors, which relate the coordinates ro on the first mirror to r on the 
second (see e.g. IVanderbei fc Traubl (|2005l )): 



or ro 



(6) 
(7) 



where c is a normalization constant chosen so ro(-R) = R — that is, both mirrors have the same 
radii. 



r (r) = ca 2 / sA m (s) 2 ds 



We also define the square root portion of Q as S(xo,yo,x,y), such that 

S(x ,y ,x,y) = V (x Q - x) 2 + (y - y) 2 + (z + h(r) - h (r )) 2 , 
Q(x ,y ,x,y) = S(x ,y ,x,y) + (h(r) - h (r )). 



(8) 
(9) 



(S will appear la t er in this derivation, and is an important factor in the S-Huygens approximation 
(jVanderbeil 120061 : iBelikov et al.l l2006al ) which is used for diffraction calculations.) The derivatives 
of Q are equal to zero w hen Q, the path length, is a constant Qq, which occurs (as shown in 



Vanderbei fc Traubl (|2005l )) when ro = ro(r), as defined by Eq. [3 and 6q = 6. 



We app roach the evaluation of Eq. [T] by using the stationary phase approximation (see ap- 
pendix III of lBorn &; Wo hi (|1999l )). which allows integrals of a certain form to be approximated for 
large values of the parameter k: 



g(x ,y )e ik ^,yo) dxQdyQ 



1 



2iria 



ky/\ a p-Y 



- T g(x c ,y c )e ik f^\ 



where (x c , y c ) are critical points where 



df_ = Q df_ 

dx ' dy 



and 



a 







d 2 f 



dx 2 
d 2 f 



d 2 f 



■E c ■ !Jc 



Xc,Vc 



a 



dx dy 

1 aP > 7 2 ,a > 0, 
-1 aP > 7 2 ,a < 0, 
i aj3 < 7 2 . 



(10) 

(11) 

(12) 
(13) 
(14) 

(15) 
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We define k "large" more explicitly as being such that fc/(xo,2/o) ^> 27r, with the assumption that 
g{xo,yo) is slowly varying and f(xo,yo) is n °t constant over the region of integration. Informally, 
the exponent must oscillate sufficiently quickly — and the amplitude g(xQ,yo) sufficiently slowly — 
that the contributions of neighboring regions tend to cancel each other out, except in the vicinity 
of certain critical points. 

By applying Eq. [10] to the equations governing the propagation through the PI A A system, we 
can determine the effect of incident wavefront profiles. (This is, of course, an approximation; the 
stationary phase approximation is the first term in an asymptotic expansion of the integral in Eq. 
[TJ The effects of higher-order terms in this expansion are discussed in Sec. 12.51 ) 

We consider first the case of an on-axis plane-wave. We can return to Eq. [TH] and rewrite its 
components as: 

2vr , . 

k = —, (16) 

E (x,y) 1 

9{xo,yo) = r = TT7T7 r for on-axis plane wave, (17) 

i\Q{x ,yo,x,y) iXQ{x ,yo,x,y) 

f(xo, 2/o ) = Q{xq, 2/0, x, y), (18) 

to make the left side identical to Eq. [TJ it remains only to evaluate the second derivatives of Q 
to determine Eq. 1121 and com plete the approximation. The algebra involved is extensive and is 
derived in full in IPueyol ( 2008 ) and Pueyo et al. ( 2011a ). and we merely report the results: for Q, 
S, ho, and h as defined above, 

a 1 

./U/3_ 7 2| ~ /, 1 ~ T ^ ' 

Vl H ' ' V ' S(x (x c ),yo(y c ),x c ,y c )*A m (ry I 

= A m (r)S(x (x c ),y (y c ),x c ,y c ), (20) 

as both and S and A m (r) are real and positive. 

The right side of Eq. [TU]thus becomes: 

!_^ =J ( I yc )e*/W = A 2mA m (r) ^O^)'^^),^,^) e 2nQ{x Q {x c ) m {y c ),x c ,y c )/X _ 
k y/\ a p - 7 2 | 6 ' 2vr m iXQ(x (x c ), y (y c ), x c , y c ) 

(21) 

If we use the fact that Q(xo(x c ), yo(y c ), x c , y c ) = Qo, we reduce this to: 

k y/\af3 - 7 2 | Qo 
In general h(r) — /io( r o) ^ S, and so: 

S(x (x c ),y (y c ),x c ,y c ) 
Qo 

e i^Q(x ,y ,x,y) dxQ dyQ ~ A m (r)e 2 ^ Q °/ A . (24) 



1, (23) 



i\Q(x ,y ,x,y) 
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and so the output wavefront for a plane wave is the desired intrinsic remapping A m (r) multiplied 
by a constant phase. As described above, A m (r) is the output amplitude distribution for ray- 
optics propagation through a PIAA system , and so, for the case of a uniform input amplitude, this 
reproduces the results of geometric optics (jVanderbei fc Traubll2005l ). 

The substitution for Eo(x,y) in Eq. [17] remains the only deviation from generality in the 
description of the input amplitude profile. We will capitalize on this to relate the amplitude 
profiles at different planes without requiring any recourse to ray optics; this is the major strength 
of the stationary phase approximation in this application. 



2.2. Case 1: known Ai(r), A m {r) 

Suppose the incident wavefront Eq(x, y) in Eq. Q]is not a plane wave, but has instead has some 
radially-symmetric apodization profile -Aj(r) with r = \/ x 2 + y 2 and zero phase: 

E ^V)= I [ ., n /' (r) re'W^'^dzodito. (25) 

J J iXQ(x ,y ,x,y) 



Then 



9(x,y) = ^, (26) 
iXQ 

and Eq. [TOl becomes 

\ 2 ™ = g(x c ,y c )e ik f^ = ^iArn^Six^MVclxc^^^^e^^ (28) 
K v / |a/3-7 2 | 2vr i\Q 

^(r^H) ^'^' 1 '^ , (29) 
and thus, given A m (r) and Ai(r), we can write 



A f (r) = A m (r)A l ^2 jT sA m (s) 2 dsj . (30) 

Note the difference between this and the product of the two remapping functions Ai(v^ and A m (v^' : 
a pre-apodization does not propagate straight through a PIAA coronagraph but is remapped to a 
certain extent. 



2.3. Case 2: known Ai(r), A f (r) 



Next, we consider the case of a known Ai{r) and Af(r). This is a design problem for a mirror 
pair: given a pre-apodizer, designed to mitigate ringing, and a final apodization profile that will 
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provide necessary contrast, what should the intrinsic profile A m {r) be? To approach this, we start 
from Eq. [30] and change variables: 



A f (r) = A m (r)Ai (c^J^ 



sA m (s) 2 dsj =A m (r)A i (r (r)), (31) 

r (r) = \l 2 j Q pA m {a) 2 da, (32) 
-A m {rf. (33) 



dr 



dr ro 

As all apodizations will be > 0, we can square both sides and change variables to convert the entire 
system to functions of r and ro: 

Mr? = ^Mro(r)) 2 , (34) 
and so ro(r) must satisfy the differential equation: 

dro rAf(r) 2 



dr r ^(r (r))2- (35) 

This is a separa b le firs t-order ordinary differential equation, whose general solution is well-known 
(see e.g. krfkenl (| 19841 )): 

/ aA f {afda= / pA0f d/3. (36) 
Jo Jo 

These cumulative functions on Ai and Af may be computed (numerically, in practice), and the 
ro(r) associated with any r determined from the Eq. [36] equality. As the apodization functions 
are both real and nonnegative, the cumulative integrals will be nondecr easing, which simplifies the 
numerics; we also recall that for a pair of mirrors with identical radius R which preserve throughput, 
f = ro(r) at r = and r = R, which should be useful as a check. A m (r) may subsequently be 
derived using Eq. [31] 

^ M = WY (37) 

We note that knowledge of A m (r) directly determines the heights of the PIAA mirrors via Eq. 
[4] and Eq. [5] and thus we can use Eq. [36] and Eq. [37] to analytically determine the mirror shapes 
required to take any amplitude profile to any other — a novel result in PIAA design. This approach 
permits a number of useful applications, such as compensating for pre-apodizer effects and chaining 
together pairs of PIAA mirrors; we return to these possibilities in Sec. [3] 



2.4. Case 3: known A m (r), Af(r) 



Finally, we consider the case of a known A m (r) and Af(r). This is a design problem for a 
pre-apodizer: given a mirror pair, what must the input apodization be to get a desired output? 
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To approach this, we return to Eq. [7J 



r (r) = cW2^ sA m {s) 2 ds, (38) 

and invert this equation (generally numerically, in practice) to find r(ro) = r^ (tq). As again A m (r) 
is real and nonnegative, the cumulative integral (and thus ro(r)) will be nondecreasing which will 
again simplify numerics. (As a check, the resulting inverse function should satisfy ^o" 1 ( r o( r )) = r 
for < r < R.) 

A change of variable in Eq. [37] then gives that 
and thus 

Mr) = (40) 

Eq. [38l and Eq. 0UJ will be valid as long as A m (r) > 0; we can easily show the mapping from 
r to tq must be invertible under this assumption by assuming the converse. Suppose there exist an 
r\ and ri such that < r\ < ri and ro(ri) = r^iri). This implies that 



cJ2 J sA m (s) 2 ds = cJ2 J sA m {s) 2 ds (41) 

=> y/^ sA m (s) 2 ds = ^ sA m {s) 2 ds. (42) 
Since A m (r) is real, and all r are nonnegative, sA m (s) 2 > and so 

rri rr2 

/ sA m (s) 2 ds = I sA m (s) 2 ds. (43) 
Jo Jo 

Since r\ < ri, 

rr2 rr\ rri 

\ sA m (s) 2 ds = / sA rn (s) 2 ds+ \ sA m (s) 2 ds (44) 

Jo Jo Jrx 

=> / sA m (s) 2 ds - / sA m (s) 2 ds = = / s^ m (s) 2 (is, (45) 

Jo Jo Jr x 

by Eq. g3j However, if A m (r) > 0, 

/ sA m {s) 2 ds> [mm A m (r)) 2 j sds (46) 

' O '1 



SdS = 2 " 2 > ° (47) 
*An(s) 2 ^ > (min An(r)) 2 (y " yj > °> ( 48 ) 



and we have a contradiction, thus the mapping will be invertible under these conditions. 
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2.5. Higher-order terms in the stationary phase approximation 

Eq. [30j Eq. [37J and Eq. HQ] are independent of any wavelengths, distances, and aperture sizes, 
and yet full diffraction modeling of a PIAA system shows the system can perform quite differently 
as these parameters are varied. Fig. [3] shows an examples of this, using mirror spacing as a variable. 
This system used both pre- and post-apodizers, which were held constant as z was changed, along 
with pair of mirrors which correspond to the same A m {r) for all values of z, and these can be seen 
in Fig. [D (This means the mirrors will always have the same shape, but scaled up and down to 
take the maximum height of the mirror from a low of 21/im at z = 1.6m to a high of 350/xm at 
z = 0.1m.) The mirror profiles were determined with Eq. l36l and Eq. [371 using the polynomial 
pre-apodizer outlined later in Sec. 13. 3[ and both a Gaussian Af(r) and post-apodizer A p (r). 



It is well-established ([Vanderbeil 120061 ) that the contrast from a pair of mirrors alone will 
be overwhelmed by diffraction from the first PIAA mirror. The inclusion of a pre-apodizer and 
post-apodizer to a mirror pair is necessary but still not sufficient to produce a system with good 
performance; diffraction modeling must be run on the actual system configuration to verify the 
system performs adequately. 

The explanation for why the system can perform so differently is that there al ways remain 



additional wavelength- and separation-dependent diffraction effects ([Vanderbeil 120061 ) which have 
not been accounted for in the design process. These residual diffraction effects correspond to 
higher-order terms in the asymptotic expansion (of which the stationary phase approximation 
is the leading order) and must be accounted for — and compensated with apodizers and system 
properties — to properly determine the expected performance of a PIAA system. A full discussion 
of these terms and their effect on PIAA systems is beyond the scope of this paper, but we can draw 
on related analysis from the laser beam shaping community to present a few salient considerations. 



The key parameter for quality of a beam-shaping optical system is /3 ([Dickey <fe Holswade 



2000 ). which in PIAA terms can be defined as 



The larger this parameter becomes, the closer a beam-shaping system (or a PIAA system) comes 
to the first term in the stationary phase approximation. For small j3 performance will be hindered; 
for sufficiently small f3, it becomes mathematically impossible to get a beam shaping system that 
works even approximate ly, as this would violate the uncertainty principle for Fourier transforms 



([Romero &: Dickeyl 119961 ) . although this will generally not be an issue for most PIAA systems. In 
practice, this means that building larger mirrors or lenses and placing them closer together will 
reduce diffraction effects; it also suggests that if (3 is sufficiently large at the longest wavelength of 
operation, the PIAA system performance will be wavelength- independent. (At the level of design, 
at least; a real system will inevitably have wavelength-dependent performance due to aberrations 
in the optics, but it should not be limited by diffraction.) 

The falloff of higher-order terms depends on the continuity class of Q(xQ,yo,x,y). For an 
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Fig. 3. — The effect of varying mirror distance on the image plane performance of a PIAA system 
with pre- and post-apodizers. Each solid curve shows the PSF of a 2-mirror s ystem with a differen t 
spacing between the mirrors, computed with the S-Huygens approximation (jBelikov et al.ll2006al ) 
to the full diffraction integral. The dashed curve is the PSF produced by the ideal Af of the system 
multiplied by a post-apodizer A p . 
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Amplitude profiles (normalized) 




0.005 0.01 0.015 

meters 

Mirror shapes 




meters 

Fig. 4. — Top. The amplitude profiles of the 2-mirror system in Fig. [3l along with the apodization 
profile of the post-apodizer and the final profile produced by the system. All are normalized to 1 
at the center. Bottom. The shape of the two mirrors, normalized to unit height on the outer edge. 
All mirror designs use an identical A m {r) and mirror radius, which gives identical mirror shapes 
scaled to greater or lesser sag following Eq. [4] and Eq. 
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infinitely differentiable Q, the higher order terms will fall off as 1//3 2 ; however, if Q has a discon- 
tinuity in its third derivative, o r Aj(r) has a discontinuity in its first derivative, then the rate of 



convergence will drop to l/yT? (|Dickev fc Holswaddl2000l ). Discontinuities in lower-order deriva- 



tives will reduce the convergence further. The input field is a plane wave with a discontinuity at 
the edge of the aperture, which immediately suggests that inserting an apodizer (see Sec. 13.2ft to 
produce a continuous Ai(r) will improve convergence. (This also suggests that the choice of A m (r) 
may affect convergence properties and thus exacerbate residual diffraction — something to keep in 
mind for design.) 

Compensating for these diffraction effects a priori is nontrivial, though a hope for future 
analysis is that the dependence of those terms on Ai(r) and A m (r) and their derivatives can be 
characterized and used to set conditions on these functions which can minimize or null the largest 
terms in the asymptotic expansion. (Ideally, this would allow PIAA systems to be designed which 
naturally minimize diffractive effects.) One fairly simple solution to mitigate some of these effects 
is to note that while the first mirror's main purpose is to redirect the light, the purpose of the 
second is to correct the phase of the wavefront so that the output has a uniform phase. The 
diffraction effects above will introduce phase errors which reduce the contrast in the final image 
plane; to compensate, we can model the residual phase at a wavelength near the center of our band 
of operation and add half this value to the surface of the second mirror. Some smoothing may 
be applied if de sired, as the residua l phase may have high-frequency components. This approach 



is also taken in iGuvon et al.l (|2010l ). to correct similar phase effects. This simple correction can 
reduce the residual light near the inner working angle by two or more orders of magnitude; the 
examples in the next section use this as a final adjustment where noted. 



3. Applications 



3.1. Post-apodizer design 



The work of IVanderbeil (|2006l ) shows that a PIAA system alone will be lim ited to ~ 10 



contra st by diffraction considerations; one method investigated in that paper and in lPluzhnik et al 



(|2006l ) to press beyond this limit was the use of apodizers preceding (a pre-apodizer) and following 
(a post-apodizer) the mirrors. The post-apodizer would do some of the beam shaping, but at the 
expense of throughput. (Pre-apodizers will be discussed in the next section.) 

The choice of the post-apodizer is a free variable; switching out post-apodizers is equivalent to 
changing the final profile of the system. The profile is selected to produce sufficient contrast at the 
image plane in conjunction with the PIAA mirrors. Eq. [30] can be modified slightly to split AAr) 
into a post-apodizer component A p (r) and the final output of the combined system A tot (r): 



A tot (r) = A p (r)A m (r)Ai 




sA m (s) 2 ds 



(50) 
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3.2. Pre-apodizer design 

One difficulty with PIAA systems is the that the beam is expected to fill the mirror; this 
produces diffraction from the edges, which (as discussed above) will limit the performance of the 
system. To compensate, we design preapodizers, which alter the discontinous step at the edge of the 
mirror into a smoother function. In general, it is better to design the pre-apodization in advance, 
since diffraction mitigation introduces an additional constraint on Ai(r) not found in A m {r) or 
Af(r); nonetheless, for completeness, we provide an example of how to retrofit a pair of mirrors 
with a pre-apodizer. 

Suppose we wish t o make a PIAA syste m which is designed to produce, as an output, a 



Gaussian apodization. (jBelikov et al.l (|2006bl ) shows these can be used to make high-contrast, 
small IWA PIAA systems.) 

A f (r) = cie _a (£) 2 ,a > 0. (51) 

Additionally, let us suppose we wish to do this with a mirror set whose intrinsic profile is a Gaussian 
which falls off more slowly: 

Amir) = c 2 e- fo (£) 2 ,0 < b < a. (52) 

The choice of a Gaussian, while providing adequate performance, is neither unique nor particu- 
larly special; we choose this particular profile for pedagogical reasons, as it produces an analytic 
representation for the preapodizer. Using Eq. we find ro(r), 



o 



C2 ^2b 



l-e- 2 K*) 

and we can write the inverse mapping for this case, as well: 



r„(r) = t/2 f ac 2 2 e~ 2h ^) 2 (53) 

(54) 



R 2 r 



^ 1 = ^o) = W^logfl-S)- (55) 



2b b V R 2 c'2 
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We can thus use Eq. 00] to write the required input apodization: 



Mr) 



C2 



-(a-6) 



Cl 
— i 

C2 
Cl Tjf 

— e 

C2 



H 2 c 



1 / 1 26r 2 \ 26 

ci A 26r 2 \^ 

C2 



i? 2 C 



(56) 
(57) 

(58) 
(59) 

(60) 
(61) 



We note, however, that most inverse mappings are not so fortunate as to have closed-form 
solutions as in Eq. [55j the inverse mapping must generally be done numerically. For example, if we 
even add a sma ll pedestal to the intr insic profile to improve the manufacturability of the mirror, 
as suggested by iPluzhnik et al.l (|2006l ) : 



(r) = c 2 



e-Hi) 2 +/ 3 
1 + /3 : 



(62) 



the mapping function becomes 
ro, P (r) = cj 



2 I P 2 .2 



4/3i? 2 



■(l-e-Ki) 2 ) 



+ 



R 2 



_ (1 + /3) 2 (l + /?) 2 6 v ; (l + /3) 2 6 

and the inverse mapping function is no longer expressible in closed form. 



(1- <.-»(*)"), 



(63) 



3.3. Mirror design 

One difficulty with pre-apodizer design is that a pre-apodizer has to meet an additional re- 
quirement that is not necessarily imposed on the other system optics: it needs to be entirely open 
across the majority of its aperture to keep throughput high, and yet has to approach zero transmis- 
sion near the edge to prevent ringing from the discontinuity in field imposed by the mirror edge. 
The pre-apodizer profile should not go all the way to zero, though, as this drives A m (r) to infinity 
in Eq. [371 hi addition, making the profile too close to or 1 at any point creates difficulties in 
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manufacturing the apodizer. (If the profile is intended to be truncated at some point, it is better 
to include this truncation in the profile Ai(r) directly so that its effects can be compensated for.) 

For these reason, it can be more advantageous to choose a desired output profile and a pre- 
apodizer which meets these requirements and design a mirror which connects the two. Consider, 
for example, a simple pre-apodizer with an open center and a polynomial falloff near the edge: 



1, r < R 



Mr) 



v 

a + 3(6-a)(^) 2 + 2(a - b) , R p < r < R (64) 

0, r > R. 



The polynomial in the region R p < r < R is chosen so that Ai(R p ) = a, Ai(R) = b, and 
dA Q^ \ r= R+ = 9A Qr^ \ r =R- = 0. Generally, we would prefer a to be as close to 1 and b as close 
to as manufacturing permit. This is not to imply that this is necessarily the best choice of pre- 
apodizer, but it has one significant computational advantage for design: the integral in Eq. I3"6l can 
be done analytically, which allows the right side of Eq. [36] to be implemented directly as a function 
instead of repeatedly determined by numerical integration. (The solution, which is a large and not 
particularly instructive 8th-order polynomial, is not shown here.) 

If we then select a final profile and a post-apodizer, Eq. [37] gives our mirror pair immediately. 
An example of this is shown in Fig. [5] using a Gaussian final profile and post-apodizer; this system 
would achieve 10 -10 contrast at 2X/D on sky. The red curve is the idealized output expected from 
geometric optics, while the blue curve includes the effects of diffraction. The parameters in this 
example are given in Tabled] the post-apodizer was proportional to e -2 ^/^ 2 , and the final profile 
was proportional to e~ 10 ( r / R ^ 2 . Propagation through the system is performed using t he S-Huygens 



approx imation to the Rayleigh-Sommerfeld integrals following the prescription in iBelikov et al 



(|2006al ). As discussed in Sec. 12.51 the shape of the second mirror is adjusted to compensate 
for phase errors introduced by the higher-order terms neglected in taking the stationary-phase 
approximation. 



Parameter 


Value 


a 


1 


b 


0.01 


R 


15 mm 


R p 


13.5mm 


z 


250mm 


A 


633nm 



Table 1: Parameters used in Fig. [5] 
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Suppression at 633nm 




5 10 15 20 25 

A/D (on sky) 



Fig. 5. — A comparison between the PSFs from a set of amplitude profiles derived from the 
equations in Section [2] ( "Perfect" ) and from the full diffractive treatment on the corresponding 
mirrors. The two-mirror system has a preselected pre- and post-apodizer, and the mirrors designed 
with Eq. [371 The second mirror is adjusted to compensate for phase errors, following § I2.51 
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3.4. Example: Gaussian input beams 

During an astronomical observation of a distant star, the incide nt field will be negligibly close to 



a uniform-amplitude plane wave. In laboratory PIAA experiments (jBelikov et al. 



2009 



Guvon et al. 



2010 



2009 



Kern et al 



Balasubramanian et al.ll2010l ). however, the starlight is simulated with a 



laser, either directly or as the output of a single-mode fiber. These beams, for lasers with high 
mode purity, will have Gaussian amplitude profiles, which effectively change the pre-apodization 
of the beam. If that beam is not sufficiently oversized with respect to the PIAA mirrors, then the 
additional preapodization from the Gaussian beam will be remapped inward and affect the inner 
working angle of the coronagraph, at least in the laboratory setting. 

Consider the case of a collimated laser beam with a single Gaussian mode, with a i-point in 
its amplitude at a radius r g : 

Mir) = c*e~te) , (65) 

entering a PIAA system with a pre-apodizer, two mirrors, and a post-apodizer. The Gaussian 
amplitude profile and the pre-apodizer combine to create an effective pre-apodization Ai(r)Ag(r), 
which is remapped to Ai(ro(r))Ag(ro(r)) by the PIAA system. Following Eq. [501 the amplitude 
profile at the output becomes 



A tot (r) = A p (r)A m (r)Ai (r (r)) A e (r Q (r)) 



A p {r)A m {r)A t icJ2 J sA m { S yds\ c e e ^ J ° . 



(66) 
(67) 



The remapped profile Ag(ro(r)) may have non-trivial effects on the inner working angle; Fig. 
[6] shows the two-mirror system from Sec. 13.31 and Fig. [5] with 4 inputs of varying beamwidth. The 
r g are shown in units of the mirror radius R, and propagation is again done using the S-Huygens 
approximation. (We note that if the remapping of the input is incorrectly ignored, the Gaussian 
profile Ai(r) will multiply the Af(r) directly; in the case of Fig. [SJ this would produce yet another 
Gaussian profile with an even faster falloff.) 

If the beam is not significantly oversized with respect to the PIAA mirrors, the effective 
change in the pre-apodization introduces significant amounts of light near the inner working angle, 
pushing it outward. This effect should be taken into consideration when preparing an input beam 
or modeling its performance in laboratory settings. 



4. Multi-PIAA systems 

One of the major problems with the use of a postapodizer is the ensuing drop in throughput, 
which removes one of the main advantages of PIAA. In addition, they sacrifice some inner working 
angle, as a post-apodizer does not get the magnification bonus of a pair of PIAA mirrors. To 
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compensate for this, we can consider using a second set of PIAA mirrors to do a remapping, with 
the intrinsic profile of the second set equivalent to the apodization profile of the postapodizer. 

To create a system with two sets of PIAA mirrors and no postapodizer, we first choose an 
overall final profile Af tSys (r) for the system, and an intrinsic A m ^{r) for the second set. We solve 
Eq. SO] to produce Ai^(r) for the second set, and as the field out of the first pair of mirrors must 
match the field into the second, we have the matching condition: 



Ai, 2 (r) = A fA (r). 



(68) 



We choose a preapodizer Aj : i(r) to minimize ringing in the system from the edges of the input 
aperture, and finally use Eq. [37] to find A mt \(r). The second mirror in the first set (and possibly 
the second set, as well) may be adjusted for phase errors. As these equations are derived from 
the leading term of the expansion of the full diffraction integral, and the appropriate choice of 
pre-apodizer and mirror size/spacing will minimize effects of higher-order terms, we can chain 
four mirrors together in a way that can provide adequate contrast even with all diffraction effects 
included. 



As an example, we can recreate the system given in 13.31 using four mirrors instead of using a 
two mirrors and a post-apodizer. The parameters chosen are identical to the ones in Table [IJ and 
the second mirror has A m ^ oc e~ 2 ( r / R ^ to match the post-apodizer. The final system profile is 
still proportional to e" 10 ^/^ 2 . The resulting intensities are shown at top in Fig. [TJ as with Fig. [5j 
the red curve is from geometric optics, while the blue curve includes diffraction using S-Huygens. 
Broadband performance is shown in Fig. [7] at bottom. 

As expected, the system gains in inner working angle, to around 1.75X/D, but its contrast is 
reduced by the additional diffraction in the system from two additional mirrors. Both this plot and 
Fig. [5] are normalized to unity at the PSF peak; the throughput of the two-mirror system is 50.9% 
on-axis, while the throughput of the four-mirror system is 87.9%. In both cases, the first mirror pair 
has sag of 87.6/im on each mirror; the second mirror pair has 86.9//m of sag. Both are shown in Fig. 
[8] (For reference, both of these designs are on flat mirror s sized in a ccordance with the low-cost /low - 
sag design principles outlined in lBalasubramanian et al.l (|2010l ) and lBalasubramanian et al.l (]201ll ).) 
We note that this approach could also be identically applied to the analysis of a pair of deformable 



(2007 


) and 


Puevo et al. 


(2011bl 



5. Summary 

In this work, we have presented three sets of equations which allow apodizations and mirror 
shapes to be designed exactly given choices for the remaining elements. We expect in the near future 
to provide explicit representations for the residual diffraction in terms of the amplitude profiles 
presented here, and the derivation of these equations from an approximation to the diffraction 
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integral provides a solid foundation for this future work. (This approximation was also shown to 
reproduce the geometric-optics amplitude distribution for a uniform input field.) A number of 
examples demonstrating the use of these equations have been provided; it is our hope that these 
will prove to be useful tools for subsequent PI A A designs. 

We have also used these tools to design and analyze systems with multiple pairs of PIAA 
mirrors, to allow significantly-improved throughput as well as smaller inner working angle. We 
expect to verify the performanc e of this class of system exper imentally in the upcoming months; 
initial progress may be found in iBalasubramanian et al.l (|201ll ). 
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Fig. 6. — The effect of an incident beam with a Gaussian amplitude profile on the system shown 
in Fig. [5j Left. PSFs from the propagation of beams with Gaussian amplitude profiles e~^'l r 9) 2 
through a PIAA system. Here R is the radius of the mirror. Right. The associated beam amplitudes, 
normalized to unit amplitude in the center. 
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Fig. 7. — Top. A comparison of PSFs between a four-mirror system with a second set of mirrors 
chosen to mimic the post-apodizer used for Fig. and the corresponding idealized amplitude 
profiles ("Perfect"). Eq. [68] is used to ensure the two pairs of mirrors couple together. The second 
mirror in the first set is again adjusted to compensate for phase errors, following Sec. 12.51 Bottom. 
The PSFs produced by the four-mirror system across a 20% bandpass; the variation in the PSF is 
quite small. 
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Fig. 8. — The radial profiles of the four mirrors in the chained PIAA system. The first mirror in 
each set is in blue, the second in green. 



